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During  the  last  six  months  of  this  contract  the  emphasis  has  been 
on  the  completion  of  studies  in  progress  for  which,  continued  support 
would  not  be  available.  The  work  has  been  concentrated  in  the  following 

three  areas:  , 

j~-  ’ , 7 7 

l)  A theoretical  and  experimental  study  of  the  O2  uptake  by  blood 
in  the  human  lung  and  in  the  artificial  membrane  oxygenator.  The  rate  of 
oxygen  uptake  by  blood  under  various  conditions  is  of  physiological  interest 
as  well  as  of  clinical  importance.  This  work  includes  a study  of  blood 
oxygenation  in  a tubular  membrane  oxygenator  under  controlled  conditions 
with  the  purpose  of  determining  the  parameters  governing  the  oxygen  transport 
under  these  conditions. 

(2)  An  improved  method  for  measuring  the  CO  diffusing  capacity  (Deo) 
has  been "developed.  Earlier  methods  havef the  shortcomings  of  requiring 
end-tidal  gas  samples  or  an  arterial  blood  sample.  A new  steady  state 
method  which  does  not  require  end-tidal  sampling  or  invasion  of  the  subject 
is  described. in  this  report 


1 


3)  Pentobarbital  is  a commonly  used  drug  as  a premedication  for 
anesthesia  and  surgery.  Hypoxia  is  always  a danger  during  anesthesia  and  4 y 't 
an  investigation  of  the  depression  of  the  ventilatory  drive  mechanisms . is- 
described. in  this  report.  It  is  anticipated  that  this  work  will  better  allow 
clinicians  to  predict  the  patients  in  whom  administration  of  this  drug  would 
result  in  depression  of  hypoxic  ventilatory  drive. 
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A.  The  Op  Diffusing  Capacity  of  Blood  at  Two  Op  Levels 

Introduction:  The  rate  of  oxygen  uptake  by  blood  in  the  human  lung  as 
well  as  in  the  artificial  membrane  oxygenator  is  of  physiological  interest 
and  of  clinical  importance.  Recently,  considerable  attention  has  been  paid 
to  the  characterization  of  blood  oxygenators  (1,3, 5, 7, 8).  In  the  present 
work,  the  study  of  blood  oxygenation  in  a tubular  membrane  oxygenator  is 
undertaken  at  various  controlled  conditions.  Since  the  diameter  of  the  tube 
is  much  greater  than  the  size  of  a red  cell,  it  is  reasonable  to  consider 
blood  as  a homogeneous  fluid  with  uniformly  distributed  red  cells  acting 
as  oxygen  sinks.  Flow  is  assumed  to  be  laminar  with  parabolic  velocity 
profile  and  the  process  to  be  diffusion  limited  with  negligible  wall  resis- 
tance. The  reaction  of  Op  with  Hb  is  considered  to  be  instantaneous;  the 
PO2  is  considered  to  be  the  same  in  the  red  cell  as  in  the  plasma.  The 
problem  is  investigated  both  theoretically  and  experimentally.  In  the  former 
part,  an  integral  method  (4)  is  introduced  to  handle  the  non-linear  diffusion 
equation.  The  main  advantage  of  the  integral  method  is  that  it  changes  the 
governing  diffusion  equation  from  a partial  differential  equation  into  an 
ordinary  differential  equation.  The  concept  of  this  method  is  to  define  a 
diffusion  boundary  layer,  6 , which  is  analogous  to  the  hydrodynamic  boundary 
layer  (6).  It  is  assumed  that  beyond  the  position  A there  is  no  mass  transfer 
and  equilibrium  exists.  By  integrating  the  governing  differential  equation 
with  respect  to  the  space  variable  over  this  layer,  a first  order  ordinary 
differential  equation  can  be  obtained  and  its  solution  is  usually  tractable. 

In  the  experimental  part,  testing  is  done  on  human  blood  at  different  oxygen 
levels  while  varying  the  rate  at  which  blood  was  pumped  through  the  system. 

By  comparing  the  theoretical  and  experimental  results,  both  the  overall 
diffusivity,  D,  and  the  diffusion  capacity,  D^,  can  be  determined.  The 
values  are  compared  with  the  resuDts  from  different  approaches  (2). 

Analysis : Fully  developed  reduced  blood  of  uniform  plasma  Op  concentration, 
C0,  enters  the  membrane  tubing  at  Z = 0 where  the  ambient  Op  pressure  is  kept 
at  constant  value  Ps.  The  steady  state  diffusion  of  O2  can  be  described  by 
the  following  differential  equation 


r75T<rl?)]  * + 


with  C = C0  at  Z = 0 

C = Cs  at  r = R,  Z > 0 

|^=0  at  r = 0,  z > 0 
dr  5 

where  D is  the  overall  diffusivity,  C represents  the  plasma  O2  concentration 
with  C = aP,  Vm  is  the  laminar  mean  velocity,  r and  Z are  the  radial  and  axial 
coordinates  Fi  represents  a sink  function  from  the  combination  of  O2  with  Hb 
with  a strength  of  Fl  = Ti(dS/dC),  where  Ti  is  the  oxygen  capacity  of  the 
normal  blood  and  S is  the  saturation. 


Equation  (l)  can  be  written  in  dimensionless  form  by  introducing  r*  = 


r*  = C ~ Co 
CS  - C0 


y = 1 - r*  yielding 
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feid-rf^i  - <l*F)y(2-y)(l-y)g- 
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with  C*  = 0 at  ? =0 

C*  = 1 at  y = 0,  f > 0 


sc 

ay 


=Oaty=l,  E>0  . 


where  Q is  the  blood  flow  rate  and 

F = 


Ti  dS 
(Cs  - C0)  dC1 


= T 


dS 


dC 


In  order  to  solve  equation  (2)  the  entire  flow  is  divided  into  two  regions 
one  being  the  entrance  region  in  which  the  diffusion  boundary  layer  develops 
and  reaches  the  center  line  of  the  tube.  This  is  followed  by  a fully  developed 
diffusion-flow  region  in  which  the  concentration  at  the  center  of  the  tube 
changes  at  different  stations  along  the  tube.  Solutions  for  each  region  have 
to  be  sought  separately. 

(A)  Entrance  region:  Let  the  concentration  profile  be  expressed  by 
a polynomial 


1 + “1  < f > - « 2<  | 


>3<  f >3 


0 < y < 6 
6 < y < l 


(3) 


where  6 is  the  dimensionless  diffusion  boundary  layer  thickness.  The 
coefficients  aq,  a2,  and  a^  are  functions  of  5 which  can  be  determined  from  the 
following  conditions 
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the  first  of  which  is  obtained  by  evaluating  equation  (2)  at  the  tube  wall, 
yielding 


-6 

ai  = FTT 


-36  2(1  + 6) 
a2  = FT T a3  = -%-+■  6 


The  only  unknown  is  6 which  can  be  solved  from  the  integral  procedure.  The 
saturation  curve  S at  pH  = 7.1*  and  37°  C can  be  approximated  by  a polynomial 


,*2 


boC*3 


S = bD  + bqC  + b2C’^  + d3i 

Integrating  equation  (2)  with  respect  to  y and  from  y = 0 to  y = 6 , and 
introducing  equation  (3)  and  equation  (4)  yields 

6?  = (1  + Tbq)  Hi(6)  + Tb2H2(6)  + Tb3H3(6)  - ^ Hl*(6) 

-?MH5(6)  - 2§goH6(6) 


00 


(5) 


3 


where  the  expression  Hi(6)  for  i = 1,  2 ...  7 are  Given  in  Appendix  A. 


Equation  (5)  gives  the  relation  between  F,  and  6 with  b]_,  b2,  b^  as 
parameters.  At  the  end  of  the  entrance  region  F = F0  where  6 has  its  maximum 
value  of  unity,  we  have 


F o = 


0.2920939  (1  + Tbq)  + O.175892  Tbp  + 0.118688  Tbq 


(6) 


Based  on  equation  (3)  the  mass  transfer  characteristics  may  now  be 
evaluated.  This  requires  consideration  of  three  distinct  kinds  of  averages: 
the  mean  value  V , bulk  mean  value  G*,  and  the  average  value  of  C*  in  each  tube, 
C*.  These  are  given  as 


= j r0*dA  = 
A J 

f VC*dA 


30  5 + 262  - 353 

10(4  + 6 ) 

56 5 - 166^  - 286 3 


JVdA 


35(4  +6") 


1126£ 


C* 


= C*d^  = (T 


6 _ riso  + 165  - 3^2  - 65 3 ~i 

io(4  + & F J 


(7) 

(8) 

(9) 


where  A is  the  cross-sectional  area  and  V is  the  laminar  velocity.  The 
saturation  at  the  exit  of  the  oxygenator  can  be  obtained  by  utilizing  equation 

(4). 

(3)  Fully  developed  region:  In  this  region  the  diffusion  boundary 
layer  loses  its  meaning  and  the  concentration  profile  is  chosen  simply  as 

c*  = fQ  - f]_y  - fay2  - f3y3  o<  y<  1 . 


By  utilizing  the  boundary  conditions 


C*  = 1 and 


3y2 


ac 

ay 


at  y = 0 


ac r 
ay 


we  have 


■ = 0 at  y = 1 

c*  = 1 - f 1 \ y + \ y2  - I y3  "j  • 


(10) 


Equation  (10)  shows  the  similarity  in  concentration  distribution  at  different 
sections  along  the  tube  in  this  region.  Integrating  equation  (2)  with  respect 
to  y from  y *=  0 to  y = 1 and  by  introducing  equations  (4)  and  (10)  yields 


where 


% = -A  In  f]_  + Bf - Efj2  + G 

A = ife  (1  + Tbl)  + + 5l4Qb^ 


„ _ 101  Tb2  . 101  Tbq 
B 720  240 


r.  _ 15997  Tbq 

b 
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The  integration  constant  G is  obtained  by  using  the  continuity  of  solution 
from  the  entrance  region  to  the  fully  developed  region,  that  is 

f l = ai  = - at  5 = 50  and  6=1 

yields  § - = _A  in  ^ + B(fx  - |)  - E(f2  - ||)  . (ll) 

Equations  (10)  and  (ll)  give  the  complete  solution  for  this  region  -where  fj 
has  its  maximum  value  of  6/5  at  5 - 50  = 0 an(i  approaches  to  zero  as 
§ - 5o  ~ °=  or  C*  — The  corresponding  expression  for  C*,  C§,  and  C*  are 

given  as 


Figure  1 presents  the  results  from  equations  (3),  (4),  (5),  (10),  and  (ll) 
which  shows  the  dependence  of  saturation  on  both  r/R  and  5.  When  5 < 1.6868 
the  concentration  profile  is  not  fully  developed,  which  corresponds  to  the 
entrance  region.  For  5 > 1.6868  the  concentration  at  the  center  of  the  tube 
starts  to  increase,  which  corresponds  to  the  solution  from  the  fully  developed 
region.  Figure  2 depicts  the  bulk  mean  saturation,  and  the  saturation 
at  the  center  line  of  the  tube  S _ . The  delay  rising  of  S is  due  to  the 
finite  interval  required  for  O2  ?o° diffuse  from  the  wall  tortRe  center  of  the 
tube.  The  dimensionles  s length  for  the  entrance  region  is  the  measure  of 
this  delay. 

Figures  3 and  4 show  the  predicted  time  course  of  blood  Pq2  at  two  ambient 
oxygen  levels.  The  pressure  ratio  Tb/Ps  is  plotted  against  the  time  ratio 
t/tog,  where  tgg  is  the  time  taken  for  T^b  to  reach  98$  P3.  It  is  clear  that 
t/tqfl  = §/?q8  irom  the  definition  of  5;  hence  the  dimensionless  time  required 
for'Tb/Ps  to  reach  98$  is  equal  to  5^0.  At  the  higher  ambient  Og  level,  a 
sharply  increasing  Fb/Ps  is  observed  in  the  early  part  of  the  process,  whereas 
Tb/Ps  increases  with  a steadily  decreasing  slope  at  the  lower  ambient  Og  level. 
This  effect  depends  on  the  range  of  the  dissociation  curve  where  the  operation 
is  taking  place.  At  the  higher  ambient  Pq2>  the  plateau  of  the  dissociation 
curve  is  reached  early  and  the  strength  of  the  sink  is  much  reduced.  After 
the  Hb  is  saturated,  additional  02  uptake  will  only  contribute  to  raising  the 
plasma  POo.  However,  at  the  lower  ambient  O2  level,  a large  sink  effect 
exists  during  the  whole  process  and  Vb/ps  then  shows  a steadily  increasing 
trend. 

The  oxygen  uptake  in  each  membrane  tubing  is 

V02  = & |T(Sb  - Sq)  + 22.4  (Cs  - C0)  Cb  } (15) 

and  the  corresponding  diffusing  capacity  becomes 
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Figure  3 

P()2  time  course  along  the  membrane  tubing  at  different  initial  saturations. 
The  value  of  tgg  is  fixed  in  this  figure;  it  is  the  time  taken  for  Fb  to 
reach ‘98%  of  Ps  when  S0  = 40%.  Holding  it  fixed  provides  a reference  value 
so  that  the  curves  for  SQ  > 40%  can  be  compared.  pH  = 7.4,  T = 37°  C. 
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Figure  4 

P02  time  course  along  the  membrane  tubing  at  different  initial 
saturations  in  a lower  ambient  02  level.  pH  » 7.4,  T = 37o  c. 
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Dl 


%2 

Ps  - ? 

Q |t  (Sb  - So)  + 22.4  (Cs  - C0)  c£  ^ 

Ps  " P 


(l6a) 


By  utilizing  the  relation  that  Q = rrDZ/2^,  then 

HDZ  |T(Sb  - S0)  + 22.4  (Cs  - CD)  c£  j- 


PL  = 


2?  (Ps  " P) 


(l6b) 


At  each  value  of  5 the  appropriate  values  of  Sb,  C*,  and  P are  all  given. 

Z is  the  length  of  the  tubing  and  D is  determined  by  comparing  the  theoretical 
and  experimental  results.  Equation  (l6a)  can  be  written  in  a different  form 
by  introducing  equation  (4)  into  it,  then 


Q STC*  ( ^ - b2C$  - 2b3C^  ) + 22.4  (Cs  - C0)  c£  j 
dGg 


tv*2 


Dl 


(16c) 


In  equation  (l6c),  dSb/dC*  represents  the  strength  of  the  sink  effect,  which 
is  inversely  proportional  to  the  value  of  Cb.  The  constant  b2  is  the  second 
derivative  of  the  dissociation  curve  evaluated  at  the  initial  condition  and 
is  always  a negative  value.  The  quantity  (dSb/det)  - b2C$)  is  by  far  the 
dominant  term  in  the  numerator.  Therefore  the  diffusing  capacity  can  be 
expressed  as 


Dloc  HU ' b2C*]  • 

The  above  expression  indicates  that  Dj,  is  mainly  proportional  to  the  blood 
flow  rate  and  the  binding  of  O2  in  the  red  cells,  which  in  turn,  can  be 
expressed  as  a function  of  the  initial  blood  condition  and  the  sink  strength 
of  the  blood  as  it  emerges  from  the  oxygenator. 

Experimental  Apparatus  and  Procedure : Figure  5 shows  the  experimental 
apparatus  which  includes  an  oxygenator,  a lucite  isothermal  testing  chamber, 
a Harvard  precision  pump  and  a 10  cc  syringe  with  water  jacket  for  blood  sample. 
The  oxygenator  consists  of  12  identical  Silastic  Medical  Grade  membrane 
tubings  (0.012  in  I.D.  x 0.025  in  O.D.  x 13*5  in  long)  with  the  entrance 
l.t  inch  coated  by  impermeable  epoxy  to  guarantee  the  full  development  of 
velocity  profile.  The  chamber  is  continuously  flushed  with  a gas  mixture  of 
desired  P02  and  4.5$  of  C02  to  maintain  constant  pH.  Heparinized  fresh  human 
blood  at  37°  C is  pumped  through  the  oxygenator  at  steady  rates  from  4.12  cc/min 
to  0.103  cc/min.  The  oxygenated  blood  sample  is  collected  and  analyzed  at 
the  exit  of  the  oxygenator  by  a SEGM  pH-Pc>2“PC02  radiometer. 

Results  and  Discussion:  Figures  6 and  7 show  the  comparison  between 
theoretical  and  experimental  results.  Each  set  of  data  is  brought  in  close 
agreement  with  the  analytical  solution  by  finding  a suitable  value  of  D.  The 
dependence  of  the  bulk  mean  saturation  Sb  on  the  initial  saturation  S0  is 


Figure  7 

Comparison  between  experimental  data  and  analytical  results  at  different 
initial  saturations  at  a lower  ambient  02  level. 
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shown  in  these  graphs.  At  two  differ' -at  ambient  pressures  Fr  = ll6  mm  Hg  and 
Ps  = 5?  mm  Hg  it  is  found  that  diffusivity  is  5.5  x 10" + cm'-ymin  respectively . 
The  difference  between  these  values  is  beyond  the  range  of  measurement  errors 
and  suggests  that  C>2  transfer  is  enhanced  at  low  O2  levels. 

Figure  8 shows  that  the  diffusing  capacity  is  a function  of  blood  fl  w 
rate  as  predicted  by  equation  (l6a) . The  higher  value  of  Dl  at  lower  ambient 
POg  and  lower  initial  saturation  is  probably  due  to  the  dependence  of  on 
the  slope  of  the  dissociation  curve,  as  predicted  by  equation  (l6c).  Because 
of  the  larger  value  of  dS/dC^  in  the  lower  part  of  the  dissociation  curve 
there  is  a strong  sink  effect  associated  with  the  lower  value  of  Ps  than 
with  the  higher  Ps.  Experimental  results  calculated  both  by  the  present 
model  analysis  as  well  as  by  the  classical  Bohr  integration  are  shown  in 
Figure  8.  Close  agreement  between  these  two  different  methods  is  probably 
due  to  the  insensitivity  of  the  pressure  gradient,  Ps  - P to  changes  in  the 
value  of  P. 

The  assumption  of  negligible  wall  resistance  in  the  present  system  is  tested 
by  oxygen  diffusion  in  plasma  where  no  sink  effect  exists.  In  this  case 
equation  (2)  is  reduced  to  a linear  differential  equation  and  the  integral 
method  of  solution  can  again  be  applied.  The  minimal  dependence  of  Op 
transport  on  y as  shown  in  Figure  9 implies  that  the  present  assumption  of 
negligible  wall  resistance  y is  justified.  The  diffusion  of  oxygen  in  blood 
whose  Hb  is  converted  to  met  Hb  is  also  studied  and  a substantial  decrease 
of  02  uptake  is  observed  as  compared  to  the  normal  blood. 

Conclusion : The  simplified  model  for  blood  oxygenation  analysis 
presented  here  is qualitatively  in  agreement  with  the  experimental  data. 

This  indicates  that  the  assumptions  of  a diffusion  limited  process  and 
parabolic  velocity  profile  are  reasonable.  Owing  to  the  flattening  of  the 
02  dissociation  curve,  the  rate  at  which  Pop  increases  along  the  capillary 
is  dependent  upon  the  ambient  O2  level.  It  is  observed  that  the  diffusing 
capacity  is  a function  of  the  blood  flow  rate  as  well  as  the  sink  strength 
and  the  initial  level  of  blood  oxygenation.  Therefore,  the  shape  of  the 
dissociation  curve  is  important  in  determining  the  value  of  D^. 

In  the  present  study  the  size  of  the  membrane  tubing  is  much  greater  than 
the  size  of  a red  cell.  For  future  work,  in  order  to  simulate  pulmonary 
capillary  blood  flow,  small  tubing  of  physiological  capillary  diameter  is 
proposed.  In  that  case,  the  parabolic  velocity  profile  is  clearly  not  valid 
and  modification  would  have  to  be  made. 


H)+(6)  = ^6  5^  + 11 6^  - 4053  + 704  ] -^(6  + 4)“'  - 8(6  + 4) 

+ 24+16  In  — ^ 4 j 


H5(6)  = -^6°  - |^65  + 7164  - + 258862  - 23056  [fi  - 8 In 


6+4 


T,  --•  --  + 4 | - 82816  i In 

H + 0 J L 


r . 4 + 6 
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B.  Measurement  of  CO  Diffusion  Capacity  Using  a New  Steady-State  Method 

Introduction;  Current  steady-state  methods  of  measuring  pulmonary  carbon 
monoxide  diffusing  capacity  while  usually  being  reasonably  accurate  and 
appreciable  have  several  disadvantages.  The  Filley  method  requires  an 
arterial  blood  sample  and  the  Bates  and  Henderson-Haggard  methods  require 
end-tidal  sampling.  The  latter  has  two  disadvantages.  In  normal  people 
shallow  and  deep  breathing  cause  falsely  low  and  high  estimates  of  Dgo 
respectively  and,  in  people  with  abnormal  lungs,  end-tidal  gas  composition 
may  not  reflect  that  of  mixed  alveolar  gas.  This  laboratory  has  recently 
developed  a steady-state  method  of  CO  diffusing  capacity  measurement  that 
does  not  require  end -tidal  sampling. 

Theory;  Symbols;  x,  xq,  xg,  X3  = inspiratory  CO,  O2,  CO2,  Infraction 

y,  Yl j Y2>  Y3  = end-tidal  CO,  O2,  CO2,  N2  fraction 

z,  zi  , Z2,  za  = mixed  expiratory  CO,  O2,  CO2,  N2  fraction 

Vp  = Dead  Space  Vol.  (BTPS) 

Vt  = Tidal  Vol.  (BTPS) 

f = frequency  (min-1)  ^ 

Dq,CO  = CO  diffusing  capacity  = D/O.863  (ml  CO-mm  Hg~ 


From  the  Bohr  relation; 

VpZ  - Vpx 
y = VT  - VD 


(1) 


The  expression  for  CO  uptake  (BTPS)  can  be  written: 


VC0  = f VT  (x  - y) 

if  R = 1 

(2) 

and 

f Vt  (x'  - z) 

if  R t 1 . 

(2a) 

Filley  et  al  have 

shown  that  (20) 

D - 

if  R = 1 

(3) 

if  R f 1. 

(3a) 

Equations  (l)  and 

(3)  and  (l)  and  (3a)  can 

be  combined  to  give; 

I = 


Vt 


vd 


if  R = 1 


where  x/(x  - z)  = I or  the  "impediment"  to  CO  gas  transfer  from  inspired  gas 
to  pulmonary  blood.  This  is  approximately  equal  to  the  reciprocal  of  the 
fraction  of  CO  removed  from  inspired  gas. 

If  R ^ 1 the  above  becomes 

x _ f Vt  + Vt  x - z 
x'-z  ~ D Vt-Vdx'-z 


I' 


(4a) 
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If  the  subject  then  breathes  at  a constant  tidal  volume  but  at  two 
different  frequencies,  the  resulting  simultaneous  equations  can  be  subtracted 
to  give : 

Iq  " 12  = (fl  " f2)  ^ if  R = 1 (5) 

with  the  last  term  "dropping  out". 

If  R / 1 and  one  assumes  that 


then 


■x.:..zciL- 

x'd)  - z(i) 


x' 


X ' Z(2) 
(2)  ' Z(2) 


0 


T . T , _ (fl  - f 2)  VT 

I1  " I2  D 


The  above  expressions  can  be  solved  for  D easily. 

As  people  exposed  to  environmental  and/or  experimental  CO  build  up  a 
significant  blood  CO  concentration  which  equilibrates  with  their  alveolar 
gas,  the  estimation  of  this  CO  "back  pressure"  becomes  important. 

If  "alveolar  back  pressure"  is  measured 

z "true"  = z "measured"  - 3/4  y (back  pressure)  (6) 

(empiric  formula) 

If  "mixed  expiratory  back  pressure"  is  measured 

z "true"  = z "measured"  - z (back  pressure)  (7) 

Thus  the  final  formula  becomes: 


dlco 


VT  (fq  - f2) 


0.863 


X 

zHtrue" 


(1) 


x 

z"true" 


(2) 


1 


Supine  subjects,  their  tidal  volume  kept  constant  by  a metering  device 
built  into  the  spirometer,  breathe  0.1$  CO  at  a high  frequency  l8  - 20/min 
for  a 4 min  rinse  period  to  achieve  a steady  state.  A 2-minute  collection 
of  mixed  expiratory  and  end-tidal  gas  is  then  made.  After  a short  interval 
a similar  rinse  and  collection  is  made  with  the  subject  breathing  at  a low 
frequency  (6-10/min)  but  the  same  tidal  volume.  All  the  above  collected 
gas  is  then  analyzed  in  a gas  chromatograph  (Beckman  GC-5,  Beckman  Instruments, 
Fullerton,  Calif.).  CO  back  pressure  is  determined  by  the  rebreathing  method 
before  the  study,  between  the  two  runs,  and  after  the  second  run  to  give  an 
estimate  of  back  pressure  for  each  run. 


Results:  To  date  three  normal  subjects  have  been  measured.  The  precision 
of  the  impediment  method  compares  favorably  with  that  of  the  Bates  method. 


abject 

DcoBabes 

% SD 

# runs 

DC0iCorr 

% SD 

# runs 

JSK 

43.3 

9.9 

8 

39.3 

7.5 

8 

SS 

29.8 

6.1 

9 

31-9 

8.0 

9 

ML 

27.G-* 

20.0 

9 

33.1 

27.9 

7 

25. 4+ 

31.3 

9 

24. C# 

32.3 

7 

*Uncorrected  for  back  pressure  fK  9 
^Uncorrected  for  back  pressure  f ^ 19 
"dJncorrected  for  back  pressure 


Discussion:  The  method  appears  to  be  theoretically  and  practically 
sound,  reasonably  accurate  and  precise,  and  it  does  not  require  end-tidal 
sampling  or  invasion  of  the  patient. 

Its  advantages  are  that  this  method  may,  in  fact,  give  the  "true  D",  a 
factor  that  may  be  especially  important  with  people  with  pulmonary  disease 
whose  "true  diffusing  capacity"  may  well  differ  significantly  from  their 
measured  diffusing  capacity.  This,  however,  remains  to  be  measured  in  such 
patients. 

Its  disadvantages  are  that:  l)  it  requires  two  separate  steady  state  runs, 
one  of  them  at  a level  of  hyperventilation  that  can  be  uncomfortable;  2)  it 
requires  CO  back  pressure  estimation.  Both  of  the  above  make  the  procedure 
time  consuming  for  patient  and  technician;  3)  it  requires  a fairly  high  degree 
of  patient  cooperation  and  therefore  might  not  be  applicable  for  ill  subjects; 
4)  it  requires  fairly  complex  equipment  i.e.  tidal  volume  regulation;  5)  the 
solution  of  the  equations  are  laborious  and  are  best  done  by  a computer;  and 
6)  it  requires  a fairly  wide  A fi  - f2>  otherwise  A Iq  - I2  becomes  a small 
number  in  which  small  errors  are  magnified  into  large  errors  in  the  final  D. 


C . Effect  of  Pentobarbital  on  Ventilatory  Driv a 

Introduction:  Although  barbiturates  are  widely  used  as  preoperative 
sedatives,  little  is  known  of  their  effects  on  the  control  of  breathing. 

While  it  is  widely  believed  that  these  agents  do  not  depress  the  ventilatory 
response  to  hypoxia  ( 15 ) , good  evidence  on  this  point  is  lacking.  In  order 
to  resolve  this  question,  we  studied  the  effect  of  pentobarbital  on  hypoxic 
ventilatory  drive  in  normal  human  subjects. 

Methods:  Seven  males  and  three  females,  all  in  good  health  between  the 
ages  of  24  and  42  years  of  age,  were  the  subjects.  All  were  born  at  low 
altitude  and  had  been  residing  in  Denver,  Colorado  (elevation  l600  meters) 
for  months  to  years.  All  studies  were  performed  in  our  laboratory  in  Denver. 
End-tidal  oxygen  and  carbon  dioxide  tensions  (Pa02  and  PpCO^)  were  continuously 
monitored.  Subjects  initially  breathed  40J&  oxygen  to  which  nitrogen  was 
added  such  that  the  P/\C>2  was  gradually  lowered  to  40  torr.  Changes  in  PaC02 
were  prevented  by  the  addition  of  carbon  dioxide  to  the  inspired  air  in 
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amounts  sufficient  to  prevent  the  development  of  hypocapnia.  Hypoxic 
ventilatory  drive  was  evaluated  as  previously  described  (see  progress  report 
No.  IV)  and  (l4)  by  fitting  the  VE  - PA02  data  points  to  the  equation: 


VE 


V0  + 


A 

PA02  - 32 


where  VE  is  minute  ventilation  in  liters  per  minute  BTPS  and  V0  is  the 
asymptote  for  ventilation  obtained  by  extrapolation.  The  parameter  A 
describes  the  shape  of  the  curve  such  that  the  greater  the  hypoxic  drive 
the  higher  the  A value.  Two  control  hypoxic  drive  studies  were  performed 
on  each  subject,  after  which  pentobarbital  2 mg/kg  was  administered  intra- 
muscularly. Measurements  of  hypoxic  drive  were  made  30,  40,  60,  75  and  in 
some  cases  105  minutes  after  injection.  Resting  ventilation  was  calculated 
as  Ve  at  PA02  = 120  mm  Hg  from  the  subjects  A value. 

Maximal  oxygen  uptake  was  measured  by  a modified  Balke  test.  With  the 
treadmill  set  at  3.5  mph,  the  subject  walked  for  10  minutes  at  a pre- 
determined grade.  Then  the  treadmill  speed  was  increased  to  7 mph  and  the 
subject  ran  at  this  predetermined  grade.  Every  2 minutes  the  treadmill 
grade  was  increased  2%.  The  test  continued  until  the  subject  signaled  his 
inability  to  continue.  Expired  gas  was  collected  in  bags  during  the  last 
30  seconds  of  each  2 minute  work  load  and  gas  samples  analyzed  by  chromatography. 

To  determine  whether  hypoxic  ventilatory  drive  was  altered  in  a given 
subject  by  pentobarbital,  the  two  control  values  on  each  subject  were 
compared  to  the  four  values  on  the  same  subject  under  the  influence  of  the 
drug  using  an  unpaired  t-test  (10).  The  significance  level  was  taken  as 
p<  .1  for  separation  of  the  two  groups.  Maximal  oxygen  uptake  in  the 
two  groups  was  compared  using  an  unpaired  t-test  (10). 

Results:  Anthropomorphic  data  relating  to  the  subjects  are  shown  in 
Table  I.  All  ten  subjects  exhibited  drowsiness  approximately  30  minutes 
after  the  pentobarbital  injection  which  gradually  dissipated  by  105-120 
minutes.  Minute  ventilation  and  end-tidal  carbon  dioxide  were  unaltered 
by  the  drug  (Table  II). 

Hypoxic  ventilatory  drive  was  depressed  by  pentobarbital  in  five  subjects 
(30-50%  depression)  and  unchanged  in  the  other  5 (Table  III  and  Figure  10). 

Within  the  group  showing  depression  of  hypoxic  drive,  changes  were  first 
observed  at  30  minutes  and  had  become  maximal  by  60  minutes.  In  all  but 
one  subject,  hypoxic  drive  had  returned  to  normal  by  120  minutes,  as  sncwn 
in  Figure  11. 

It  appeared  that  this  group  of  subjects  might  be  comprised  of  2 separate 
subgroups  with  different  responses  to  pentobarbital.  In  these  2 subgroups 
ages  were  not  statistically  different  nor  were  weights,  resting  ventilation, 
or  PAC02.  There  were  3 females  and  2 males  in  the  depressed  group  and  5 
males  in  the  nondepressed  group.  Resting  heart  rates  were  lower  (see  Table  III) 
58  + 3 as  opposed  to  71  + 5 beats  per  minute  (P<.05)  in  the  group  that 
showed  no  change  in  hypoxic  ventilatory  drive.  In  addition,  maximal  oxygen 
uptake  was  higher  in  this  group,  (49  + 2.63  as  opposed  to  33  + 2.05)  (p<.00l), 
(Table  IV) . 
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TABLE  I.  ANTHROPOMORPHIC  DATA 


• 

Sex 

Age  (yrs.  ) 

Ht.  (cm) 

V/t.  (kfi) 

CAH 

F 

29 

155 

53.  6 

REM 

M 

31 

189 

64.  7 

AT 

M 

26 

174 

70.  0 

ML 

F 

24 

160 

54.4 

IIC 

M 

37 

170 

57.  1 

JSK 

M 

24 

180 

58.  9 

i'H 

M 

28 

175 

78.  0 

FIH 

M 

33 

188 

97.4 

LL 

M 

41 

175 

82.  5 

BK 

F 

42 

165 

70.  0 

Mean 

31.  5 

173.  1 

68.  7 

SEM 

2.  09 

3.  50 

4.  45 

25.77  .85  .53  .92  IP. 68  .615  .45  .90  13.92 


TABLE  II.  Ccn 


•79  .63  1.05  11.86  .63  ,9S  2.06  42.92 


ABL6  1 1 T • HYPOXIC  DUTVK  DURING  PENTOPAKBITAL  TN  2 SUBGROUPS: 


depression 


Subject 

Posting 

Heart  Kate 
He.itr/min 

V°2  Max 
r1  't:g/ni  n 

0 

Time  in  min. 
30  4 5 

after  Pcntobarbi tal 

r.o  7 5 

1 20 

CAII 

80 

32.3 

358  331 

216 

149 

162 

168 

333 

ML 

63 

34.6 

135  151 

154 

154 

98 

08 

126 

1IC 

76 

37.2 

81  85 

48 

30 

63 

73 

123 

LL 

00 

35.6 

199  250 

131 

144 

3 

79 

111 

BK 

55 

25.5 

100  97 

69 

62 

44 

56 

99 

Mean 

SEM 

70. 8+ 

5.03 

33. 4 + 
2.05 

170.8 

48.21 

123.6 

30.18 

103.8 

22.62 

74.0 

26.82 

92.8 

17.63 

167.2 

42.9 

r<.os 

+ P<.01+ 

• 

110M 

DEPRESSION 

Resting 

Heart  Rate 
Subject Beats/min 

V0-,  Max 

r.l/kn/nin 

0 

Hypoxic  Ventilatory  Drive  (A) 
Time  in  min.  after  Pentobarbit 
30  45  60 

OS 

al 

75 

REM 

52  . 

55.0 

149  121 

150 

105 

132 

119 

AT 

52 

40.8 

157  166 

128 

179 

141 

119 

SK 

54 

54.0 

175  130 

223 

209 

166  ' 

1:1 

FHu 

68 

40.  3 

53  56 

68 

55 

60 

55 

FHei 

63 

47.6 

166  133 

121 

166 

14  3 

144 

Mean 

SEM 

57. 8+ 

3.26 

49.  r 
2.63 

131.4 

19.59 

138.0 

25.16 

150.8 

26.87 

128.4 

18.00 

115.6 

16.04 

Paired  T 

test  was  used 

to  conpare  2 

control 

A values  v/ith  4 A values  after 

rent.oharbital 

r 


k'5 


TABLE  IV 

Maximum  O2  Uptake  ml/kg/min 


Depre 

used  Group 

Non-Depressed  Group 

Subject 

VO?  max. 

Subject 

VOp  max. 

SK 

54.0 

CAH 

32.3 

REM 

55.0 

ML 

34.6 

AT 

48.8 

HC 

37.2 

F Hei 

47.6 

BK 

25.5 

F Hu 

40.3 

LL 

35.6 

X 

49.1  + 2.63 

X 

33.4  + 2.05 

p<  .001 


EFFECT  OF  PENTOBARBITAL  ON  HYPOXIC 

VENTILATORY  DRIVE 


HYPOXIC  VENTILATORY  DRIVE  DURING  PENTOBARBITAL 
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Discussion:  In  healthy  subjects  pentobarbital  in  amounts  commonly  used 
to  produce  mild  sedation  has  no  effect  on  resting  ventilation  or  resting 
PAC02  but  did  depress  hypoxic  ventilatory  drive  in  half  the  subjects. 

Absence  of  alterations  in  resting  ventilation  is  in  agreement  with  Stephen 
et  al  (13)  who  found  no  hypoventilation  with  a comparable  dose  of  secobarbital. 
Keats  et  al  (12)  found  an  increase  in  ventilation  after  1.4  mg/kg  pento- 
barbital but  a decrease  after  2,8  mg/kg.  The  most  clear-cut  evidence  of  a 
respiratory  effect  of  pentobarbital  was  a depression  of  hypoxic  drive  in 
five  of  our  subjects  and  no  change  in  hj'-poxic  drive  in  the  remaining  five. 
Little  has  been  written  regarding  the  effects  of  barbiturates  on  hypoxic 
ventilatory  drive.  In  fact,  the  only  previous  study  found  (ll)  demonstrated 
in  one  subject  that  triple  conventional  doses  of  cyclobarbital  resulted  in 
some  depression  of  ventilation  during  hypoxia.  However,  only  very  mild 
hypoxia  was  employed  (P^Op  = 60  mm  Hg) , which  is  insufficient  to  produce  a 
well  developed  hypoxic  response  in  normal  man  (l4) . 

In  half  the  subjects,  there  was  no  depression  of  hypoxic  ventilatory 
drive  with  pentobarbital.  There  appeared  to  be  two  populations  represented 
in  our  sample  which  differ  with  respect  to  the  suppressibility  of  their 
hypoxic  drive  with  pentobarbital.  These  subgroups  also  differ  in  their 
physical  condition  as  measured  by  maximal  oxygen  uptake  and  resting  heart 
rate.  Subjects  whose  ventilatory  drive  was  unaffected  by  pentobarbital  and 
who  were  all  involved  in  athletic  activities  on  a regular  basis  had  lower 
resting  heart  rates  and  higher  maximal  oxygen  uptake  than  the  subjects  in 
the  group  that  showed  depression  of  hypoxic  drive.  Byrne-Quinn  et  al  (9) 
found  that  hypoxic  drive  was  lower  in  true  athletes  at  rest  than  in  nonathletes. 
While  we  found  no  significant  differences  in  control  levels  of  hypoxic 
ventilatory  drive  between  the  two  groups,  some  other  factor  distinguishing 
athletes  from  nonathletes  could  account  for  the  differences  in  response  to 
the  drug.  Body  weight  was  similar  in  the  two  groups  and  moreover  the  drug 
was  administered  on  a weight  basis. 

However,  it  is  possible  that  differences  in  lean  body  mass  or  in  drug 
metabolism  or  redistribution  could  account  for  this  difference.  Unfortunately, 
we  do  not  have  barbiturate  blood  levels  on  any  of  our  subjects  and  cannot 
exclude  the  possibility  that  differences  in  drug  metabolism  may  have 
contributed  to  the  different  effects  in  the  two  groups.  This  study  was 
not  carried  out  with  a double  blind  design  because  we  believed  that  our 
subjects  would  all  have  been  able  to  readily  distinguish  the  difference 
between  pentobarbital  and  placebo. 

This  study  suggests  that  sedentary,  less  physically  fit  individuals 
are  more  likely  to  exhibit  depression  of  hypoxic  ventilatory  drive  where 
pentobarbital  is  administered. 
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